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ABSTRACT 

At the LHC, where hard multi-jet events are common, large imbalances in the 
total jet-p T can arise due to the interaction of jets with the plasma formed in heavy-ion 
collisions. We find that this missing-p T spectrum yields a measurement of the initial 
temperature of the plasma. Using commonly accepted guesses for this temperature, and 
ISAJET estimates for jet cross sections, we predict an order of thousand events at LHC 
with unbalanced p T > 50 GeV. 



TIFR/TH/94-29 
February 2008 



For the heavy-ion collider mode of the LHC, with \fS — 5.5 TeV, the cross section 
for multi-jet events in which none of the jets is less energetic than, say, 50 GeV is of the 
order of several //b. In a year's run one expects to see about 10 6 -10 7 jet events. Of these 
about 10 4 events contain four or more jets. The expected distribution in the number 
of jets, N, is shown in Fig. 1. Since the underlying event comes with a multiplicity 
density, dN^/dy, of 2000-3000 pions per unit of rapidity, the transverse energy within 
a cone of AR = 0.5 is about 20-30 GeV (assuming (p T ) 0.5 GeV per pion). Thus 50 
GeV jets are easily observed above this background. 



1000 r 



100 - 



Figure 1 

The distribution in the number of jets obtained from a sample of 10 5 events gener- 
ated in ISAJET with the kinematics mentioned in the text. 

If simultaneously, a quark-gluon plasma is formed, then the evolution of the jet 
is coupled to the plasma through strong interactions. This opens up the possibility of 
measuring plasma properties using jets. Such measurements could truly be called hard 
probes of the plasma. The simplest measurements are kinematic, involving changes in 
the 4-momentum of jets as they evolve through the plasma. It should be emphasised 
that these are strictly probes and not signals of the plasma. 

There is a history of related studies for dijet events. The first suggestion [1] was 
to study acoplanarity in 2-jet events. A static model of the plasma was used to show 
that small acoplanarities can exist. This was extended in [2] to include the effects of 
hydrodynamics and a mixed phase. In [3] it was shown that a hadron gas could give 
numerically similiar results. Acoplanarity is an useful variable only for 2-jet events. 
Very recently energy-energy correlations in 2-jets events have also been investigated [4] . 
A parallel line of attack is to study the energy loss of a test parton evolving through 
the plasma. This gives rise to the observation of jet quenching [5] . 

In the context of forthcoming LHC experiments, we believe that the simplest mea- 
surement is of missing transverse momentum, j> T . This is a standard tool for many 
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other kinds of physics, and one of the LHC detectors can easily perform this measure- 
ment. Such a measurement need not be restricted to dijet events. If an experiment 
with 4ir coverage observes exactly AT-jets, with a given $ T , then this measurement can 
be converted into a statement about the initial temperature, T , of the plasma. The 
experimental analysis is made robust by three pleasant aspects of the signal. 

i) The j4 T -spectrum is almost independent of the initial jet energies when each jet is 
harder than 50 GeV. Thus there is no need to guess the energy of the jet at the 
hard vertex, before its interaction with the plasma. 

ii) The j4 T -spectra for different iV are related. Taken together, they provide a cross 
check for the physical origin of the missing momentum. Unlike acoplanarity, this 
measurement is not restricted to the dijet sample. 

iii) Using well-accepted guesses for T at LHC, we estimate about a thousand events 
with j) T > 50 GeV. This would yield an accurate experimental measurement of T . 

Backgrounds coming from lack of full 4tc coverage, contamination by the underlying 
event, and the production of neutrinos and non-standard- model particles can be made 
small, and will be discussed at the end of this letter. 

In this letter we outline the computation of the ^ T -spectrum, and show the relation 
between the ^ T -spectra for different number of jets. After this we proceed to display 
our results. We find that the ^ T -spectrum is insensitive to the initial energy in the jet, 
but strongly sensitive to To. It is to be emphasised that j) T arises when considering only 
the jets. If the p T of all hadrons in each collision were to be added up, there would be 
no j) T . 

The Computation: The partons which give rise to jets are produced through a 
hard scattering in the earliest times of heavy- ion collisions. In the absence of a medium, 
the measurement of the momenta of the particles constituting the jet would give the 
momenta of these partons. This enables us to compute jet production cross sections 
in perturbative QCD. In the presence of a medium, final state scattering would occur. 
If the medium is a quark gluon plasma, then the individual scatterings may again be 
computed in perturbative QCD. 

The energy of the constituents of the plasma are typically of the order of a GeV, 
whereas we shall consider jets carrying energy of order 50 GeV. The difference in the 
scales translates to the kinematic statement that the fractional momentum change of 
the jet in each scattering is small. Hence, the individual scatterings may be considered 
independent. We can then write the ]6 T -spectrum in an A^-jet process as the convolution 
of many re-scattering probabilities, Fj, — 



V N (i> T ) 
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Here the cross section a may be considered to be differential in as many other variables 
as required. The individual distributions, F i: may also be considered to be functions of 
these other variables. Note that the distribution in j4 T is solely given in terms of the 
interaction between the jet and the plasma. The details of the hard scattering process, 
including the perturbative expansion of the matrix elements, possible resummations of 
initial state radiation, jet definitions, etc., are all subsumed into the factors of a. 

The probability of generating a missing transverse momentum fa T for a hard parton 
of species j is given by 



where the index i runs over the species of partons in the plasma and Ni{r) is the 
number density of partons of species i at time r. The hard parton is assumed to 
interact with the plasma for times 7~o < r < Tf. The initial time, To, will be taken 
as the thermalisation time of the plasma. The final time, Tf, can be obtained from 
kinematic considerations. The function contains the singular part of the invariant 
cross section for the scattering of partons of species i and j. Since the relevant jets are 
gluonic, in the rest of this paper we shall lighten the notation by dropping the index j. 

A two-dimensional Fourier transform decouples the 5- function in Eq. (1) and yields 
the characteristic function of the ^-distribution as the exponential of the characteristic 
functions of the distribution F [6]. Thus, after normalisation, we obtain 



We have used the fact that /j is a function only of the magnitude of j4t to perform 
the angular integration in the Fourier transform and obtain the Bessel function Jo(b$ T ). 
The subtraction in Eq. (4) makes F(0) = 0, i.e., S(0) = 1, and hence takes care of the 
normalisation of V N (fi T ). The lower limit, /j,, of the integral over q regulates a leftover 
logarithmic divergence, and is chosen to be the Debye screening mass evaluated at one- 
loop order. The ^-distribution is recovered by an inverse Fourier transform, which, 
after the angular integration is performed, takes the form 





S N (b) = / d 2 i> T exp[-ib-i> T ]V N (t> T ) = exp[F N (b)}, 



(3) 



where 




(4) 




(5) 
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Note that many variables are implicit in this formula. The dynamics of the plasma 
enters through the distributions Mi and explicitly also through tq. The kinematics of 
the jets enter through the quantities E a and tJ. We specify these dependences next. 

We begin with the evolution of the plasma. In this letter we assume the plasma 
to be described by longitudinal boost-invariant hydrodynamics. We shall take as our 
initial conditions a head-on collision of two nuclei of mass number A. The plasma 
volume is then cylindrically symmetric, with radius R = roA 1 / 3 (ro = 6 GeV -1 ). A 
local temperature may be defined, T(r). The particle distributions, as functions of this 
temperature, have an expansion in a s . To leading order, they are 

Mi(r) = ^T 3 (r), where 9i = 4 _£^ 9 _ (6) 

The evolution of the local temperature keeps the entropy density fixed. Hence 

tT 3 (t) = r T 3 . (7) 

Tq may be related to the pion multiplicity density assuming that there is no entropy 
generation during hadronisation of the plasma. We can write 

ToT ° - 3C(W^' { ) 

In most of our calculations, we shall use this formula to specify T in terms of the 
multiplicity. However, in actual experiments, the jet measurement would, among other 
things, check the validity of this formula. 

The functions /j also have expansions in a s . Again we retain the leading term, and 
obtain 

m ^ u. _ SI, tort = g 



MQt) = hi^^, where K = *> ^ ." y _ (9) 
Qt I 2, for % = q,q 

From Eqs. (6) and (9) we obtain 

F (9riT ) = A%1, where A= 36(4 + y C(3) T'(r). (10) 

We have used the one-loop order expression for a s (q T ). 

Next, we define the kinematics of the multijet event. Assume the hard vertex to 
be at distance r from the axis of symmetry. The momentum of the a-th parton can be 
taken to be 

p a = £ Q (l,sechy Q cosr,sechy a sin6> a ,tanhy Q ), (11) 
A simple computation now yields 

T f = ^r 2 + R 2 -2rRcos<j), where sin(0 - 0° ) = ~ sm6 a . (12) 
J R 
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Since the polar angles and the position of the hard vertex, r, are not measured, 
these variables must be averaged over. This is the final ingredient in the computation. 
We assume that the cross section for a hard process is just the convolution of the hard 
cross section with a nuclear density function. Assuming an uniform density, we obtain 

v N (M = ±J dbbJ (t T b) J drr ^ R3 r2dn s N (b). (13) 

The volume element (dfl/Q) averages over the jet angles 9 a (a = 1, • • • , N) with uniform 
measure. The distribution is clearly a function of the energies and rapidities of the jets. 
The ^ T -spectrum is now completely specified. 

Some features: A few remarks may help understand the results of the numerical 
work reported later. In this approach to the interaction between jets and a plasma, each 
jet undergoes successive scatterings from the soft plasma independently of the others. 
The net effect is a large imbalance in p T . This remaining p T is found in the soft hadronic 
background underlying the jets, but is invisible in the jet measurement. 

Although each single collision of the jet with the plasma only changes the vectorial 
p T by a number of the order of T, the number of collisions along the path is large. It is 
proportional to Af integrated along the path, and hence increases at Tq . Note that (jfe) 
computed from Vi($ T ) equals that from F(fi T ). Resummation only changes the shape 
of the ^ T distribution, introducing a longer tail at the cost of the small-^ T region. 

In [1] QCD initial state radiation was also resummed a la Collins and Soper [7] 
and added to the ^-distribution. We have not done this for a simple reason. Such 
a resummation is performed for an A-jet inclusive cross section. For large j4 T this 
estimates the effect of events where the number of hard jets is larger than A. Since we 
consider exclusive A-jet events, and large tf T , inclusion of such a term in our formulae 
would be erroneous. 

In [2] a component of the scattering from the mixed phase was also included. It 
was observed in [4] that jets escape the plasma before the mixed phase can form. The 
kinematics in Eq. (12) also lead to the same conclusion. Hence we have not included 
mixed phase cross sections in our computations. A further ingredient in some studies 
of the plasma is the inclusion of transverse hydrodynamics. While this decreases the 
plasma lifetime, it also enhances the plasma volume. Thus the net effect on the f> T - 
spectrum may be small. It will be studied in more detail in a later paper. 

Results: For orientation we display results for a single jet travelling through a 
plasma. This may arise, for example, as the balancing jet for high-p T Drell-Yan events. 
In this case f T is the imbalance in p T between the jet and the lepton pair. For A = 1 
dQ/Q = (16/277. We take A = 200, to = 0.1 and the jet rapidity to be zero. For 
dN^/dy = 2000, the ^ T -spectrum is shown in Fig. 2. We obtain (^|) = 16 GeV 2 , and 
find that the distribution goes to zero rapidly near j) T = 40 GeV. The spectra for jet 
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energies in the range 100 to 1000 GeV are almost equal. The dependence on A at fixed 
T is small; (jf>%) changes by 6% in going from A = 50 to A = 200. Since this change 
is due only to the thickness of the plasma firetube, it also demonstrates that we can 
neglect changes in ^-distributions for small changes in the impact parameter of the 
colliding nuclei. Increasing dN n /dy to 3000 broadens the spectrum a little; (jf>%) = 18 
GeV 2 in this case. However, it is unlikely that such small j> T can be reliably measured 
above the soft hadronic background. 
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Figure 2 

Vi(j T ) plotted against $ T (in GeV) for dN^/dy = 2000. 
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Figure 3 

P>(p c T ut ) plotted against p£ ut (in GeV) for multiplicity densities of (a) 2000 and (b) 
3000, per unit of rapidity. The curves are for 2 (full line), 3 (dashed line) and 4 (dotted 
line) jets at zero rapidity with = 400 GeV. 



Multijet events display a much harder ^ T -spectrum. At fixed dN^/dy = 2000, the 
values of {f T ) are 23 GeV 2 , 29 GeV 2 and 33 GeV 2 for iV = 2, 3 and 4, respectively. In 
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Fig. 3 we display 

„cut 
fPT 

r>(p c ^) = i- df T v N (i> T ), (w) 

Jo 

i.e., the fraction of events with j4 T > p^ ut . For dN^/dy = 2000, we see that the fraction 
of events surviving a cut of j4 T > 50 GeV are 1%, 5% and 10% of the total. For 
comparison, the corresponding numbers increase to 3%, 10% and 17% respectively, 
when dN^/dy = 3000. Thus, the distinction between different values of T is easily 
observable, and becomes easier with increasing N. 

Backgrounds: The main background is the uncertainty in j4 T measurements due 
to the high multiplicity in the underlying event. If the transverse energy within a jet 
cone of AR = 0.5 due to this effect is E®, then the ^ T uncertainty in an iV-jet event 
becomes E°VN. Since E° is 20 GeV (30 GeV) for dN^/dy = 2000 (3000), the ^ T 
uncertainty ranges from 28 GeV (42 GeV) for iV = 2 to 40 GeV (60 GeV) for N = 3. 
This background can easily be eliminated by observing distributions with a cut j4 T > 50 
GeV. As noted in the preceding paragraph, the number of events is still very large. 

A second important background comes from limited angular coverage of the detec- 
tor. Missing p T could be generated by a jet lying outside the detector. An ISAJET 
simulation showed that for a detector coverage of \y\ < 5 and all polar angles, the back- 
ground to signal ratio can be reduced to less than 10 -5 by requiring that the observed 
jets all lie within \y\ < 2. 

A third background comes from W or Z production with a subsequent leptonic 
decay giving a neutrino which escapes undetected. This background is only 4% of the 
signal, since the cross section is damped by the ratio of the weak and strong couplings 
times the semi-leptonic branching ratio of the vector boson. It may be reduced even 
further by vetoing on a hard lepton in the direction of ]6 T . Kinematically similiar 
backgrounds may arise from production of non-standard-model particles. These are 
small, and one could explore ways of reducing them through topological or kinematic 
cuts. All these backgrounds will be studied in greater detail in a separate publication. 

The suggested experiment: For an LHC detector with full polar angle and 
\y\ < 5 coverage, we suggest the following measurement: 

1) Select jet events with all jets harder than 50 GeV and \y\ < 2, and bin according 
to the number of jets N. 

2) Veto events with missing transverse momentum less than 50 GeV. 

3) Veto events with one or more hard charged leptons in a cone AR — 0.5 around the 
direction of the missing momentum. 

4) Match the j> T distribution against the formulae in this letter, varying Tq for each 
N. The result is a measurement of To. 

About 10 5 events are expected to pass the cuts in steps 1-3 for an integrated luminosity 
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of 10 4 ^b" 1 . Typically 10% of the events will have j> T > 50 GeV. If heavy -ion collisions 
give rise to a quark-gluon plasma at the LHC, then in a year of running one could expect 
to obtain a good measurement of the initial temperature of the plasma. 
I would like to thank Nirmalya Parua for help with ISAJET runs. 
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